On the full quantum trispectrum in multi-field DBI innation 
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We compute the leading order connected four-point function of the primordial curvature pertur- 
bation coming from the four-point function of the fields in multi-field DBI inflation models. We 
confirm that the consistency relations in the squeezed limit and in the counter-collinear limit hold as 
in single field models thanks to special properties of the DBI action. We also study the momentum 
dependence of the trispectra coming from the adiabatic, mixed and purely entropic contributions 
separately and we find that they have different momentum dependence. This means that if the 
amount of the transfer from the entropy perturbations to the curvature perturbation is significantly 
large, the trispectrum can distinguish multi-field DBI inflation models from single field DBI inflation 
models. A large amount of transfer Tns S> 1 suppresses the tensor to scalar ratio r oc T^J and the 
amplitude of the bispectrum oc T^J and so it can ease the severe observational constraints on 

the DBI inflation model based on string theory. On the other hand, it enhances the amplitude of 
the trispectrum r^" oc T^s/nx" 2 for a given amplitude of the bispectrum. 
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Precise measurements of the cosmic microwave background (CMB) anisotropics such as those obtained by the 
WMAP satellite [l[ provide valuable information on the very early universe. Any theoretical model that attempts to 
explain the evolution of the universe before the big bang nucleosynthesis will also have to explain the observed CMB 
anisotropies. Even though these anisotropies are almost Gaussian, a small amount of non-Gaussianity is still allowed 
by the data @, S m an d future experiments such as PLANCK [5] might detect this small amount of the primordial 
non-Gaussianity. Even though the simplest slow-roll single field inflation models predict that the non-Gaussianity of 
the fluctuations should be very difficult to be detected recently, theoretical models which can produce sizeable 
non-Gaussianity has been extensively studied by many authors 
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Among them, the Dirac-Born-Infcld (DBI) inflation model can produce large non-Gaussianity by the fact that the 
sound speed of the perturbations can be much smaller than one due to the non-trivial form of the kinetic term [67t |68| . 
DBI inflation models are also well motivated by string theory [ID, H3, [til, [zO, [zE [ll] . In this model, the inflaton fields 
are identified with the positions of a moving D3 branc in the six-dimensional internal space. The dynamics of the D3 
brane is described by the DBI action. However, recently it has been pointed out that single field DBI inflation driven 
by a mobile D3 brane with large non-Gaussianity might contradict with the current WMAP data. For current and 
stringent observational constraints and consequences of DBI-inflation see (73l . l74l . [tBL l76l \tA l78l l79l l80l . l8ll . [H, HU . 

One way to avoid these constraints is to consider multi- field DBI models [36|. Since the position of the brane in 
each compact direction is described by a scalar field, DBI inflation is naturally a multi-field inflationary model [84j |. 
As first pointed out by [8o| . in multi-field inflation models, the curvature perturbation is modified on large scales 
due to the entropy perturbation. Even though there are some works considering multi-field inflationary models with 
non-canonical kinetic terms [45|, H(| H3, HI], the consistent analysis for the entropy modes in the multi- field DBI 
inflation model has started only very recently [H, [36|, H3] • 

In [H,[36;, 37], the three-point functions in the small sound speed limit and at leading order in the slow-roll expansion 
were obtained and it was shown that in addition to the purely adiabatic three-point function, there exists a mixed 
component (Qa-(k.i)Q s 0^2)Qs0^3)) where Q a and Q s are the adiabatic and the entropy perturbations, respectively. 
Since the momentum dependence of the three-point function from the adiabatic modes was shown to be identical 
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with the mixed component, the shape of the bispectrum of the curvature perturbations remains the same as in the 
single-field case, while the amplitude is affected by the entropy perturbation. 

Even though previous works on the non-Gaussianity in the multi-field DBI inflation model are mainly limited to the 
bispectrum, it is expected that the cosmic microwave background (CMB) trispectrum also provides strong constraints 
on early universe models. At the moment, the constraints are rather weak given by |tjvx| < 10 s [89. |90|. where t^l 
denotes the size of the trispectrum. However, PLANCK will tighten the constraints significantly reaching \tnl \ ~ 560 
[qTI | . Although these estimations are obtained assuming local type non-Gaussianity, we expect similar constraints can 
be obtained for non-Gaussianity from DBI inflation. It is also worth noting that the analysis in the single field DBI 
inflation model shows that the trispectrum from the contact interaction diagram is enhanced in the small sound speed 
limit as t^l ~ l/ c s [HI]. Even though the contribution from the scalar exchange interaction was overlooked in this 
analysis, recently, it has been confirmed that the same scaling holds for the trispectrum from this interaction [62l[63|. 
As in the bispectrum case [Hj], the observational constraints depend on the shape of the wave vectors' configuration. 
Therefore, it is important to calculate the full shape dependence of the trispectrum in the multi-field DBI inflation 
model. For the details of the observations of the CMB trispectrum, see [§3, 0, HE] • 

For this purpose, recently we have calculated the quantum trispectrum from the contact interaction and obtained 
the simple understanding of the origin of the interaction terms of the multi- field DBI inflation [64j . In this paper, as a 
natural continuation, we will calculate the scalar exchange trispectrum and obtain the complete theoretical prediction 
for the quantum trispectrum of multi-field DBI inflation. For a related work, where only the purely entropic component 
of the quantum scalar exchange trispectrum was calculated, see version three of .55]. 

The structure of this paper is as follows. In the next section, we shall introduce the model and some basic notation. 
In section HIT] we summarise the results obtained in [64| for the contact interaction trispectrum. After calculating the 
scalar exchange trispectrum in section ITVl the momentum dependence of the trispectra coming from the adiabatic, 
mixed, and purely entropic contributions are studied separately in section |Vj The observational constraints on the 
multi- field DBI inflation model are discussed in section PVTI Section IVlII is devoted to the conclusions. 



II. THE MODEL 



In this section, we will introduce the multi-field Dirac-Born-Infeld (DBI) inflationary model. We will present the 
background evolution equations and define some basic notation. 

The multi-field DBI inflation model is described by the following action [96[ 



S 



d^x. 



R 



f 2P(X,^) 
^l-2f(4> I )X-l)-V(ct> 1 ) 



(1) 



where we have set 8ttG = 1, R is the Ricci scalar, (f) 1 are the scalar fields (/ = 1, 2, iV^), /((f) 1 ) and V(4> 1 ) are 
functions of the scalar fields determined by string theory configurations and X is defined in terms of the determinant 
V = det(<5£ + fGud^tfdu^) as X = (1 - V)/(2f). Here G u is the metric in the field space. We assume that P 
is a well behaved function of (f) 1 and X. It is also shown that X is related to the kinetic term of the scalar fields as 

dm 



X = G U X IJ - 2fx} I Xf ] 



ifxfx/Xp - 8/ 3 x/A/A/A L L 



X 



IJ _ 



Xj — GikX 



KJ 



(2) 
(3) 



where the brackets denote antisymmetrization. It is worth noting that even though X = X ( = GijX ij ) in the 
homogeneous background, this does not hold if we take into account the inhomogeneous components. 

In the background, we are interested in flat, homogeneous and isotropic Friedman-Robertson- Walker (FRW) uni- 
verses described by the line element 



ds 2 = -dt 2 + a 2 (t)Si j dx i dx j , 
where a(t) is the scale factor. The Friedman equation and the continuity equation read 

3H 2 = E , 



(4) 



(5) 



E = -3H (E + P„) , 



(6) 
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where the Hubble rate is H = a/a, a dot denotes derivative with respect to cosmic time t, Eq is the total energy of 
the fields which is given by 

E = 2Xi J P , XIJ - Po, (7) 

and the subscript denotes that the quantity is evaluated in the background. 

For this model the speed of propagation of scalar perturbations ("speed of sound"), c s , is given by 

c2=(, . (8) 

Since we are interested in the inflationary background, we assume that the form of /((j) 1 ) and V{4> 1 ) are chosen 
so that inflation is realized at least for 60 e-foldings. In order to characterize this background, we define the slow 
variation parameters, analogues of the slow-roll parameters, as: 

We should note that these slow variation parameters are more general than the usual slow-roll parameters and that 
the smallness of these parameters does not imply that the field is rolling slowly. We assume that the rate of change 
of the speed of sound is small (as described by x) but c s is otherwise free to change between zero and one. 

We shall consider perturbations on this FRW background. We decompose the scalar field cj) 1 into the background 
value 4>q and perturbation Q 1 in the flat gauge as, 

<f> I (x,t)=4, I (t)+Q I (x,t). (10) 

Furthermore, as was done in [93], we decompose the perturbations into instantaneous adiabatic and entropy per- 
turbations, where the adiabatic direction corresponds to the direction of the background fields' evolution while the 
entropy directions are orthogonal to this. We introduce an orthogonal basis e l n , with n = 1,2, N^, in the field space 
so that the orthonormality condition are given by [35j ] 

1 1 - c 2 

e n e m i — — S mn -SmiSni , (11) 

Cs C s 

where the adiabatic basis is defined as 



(12) 

Notice that the length of the basis vector is c s and that of the other basis vectors is l/c s . If we consider the 
two-field case (I = 1,2), the field perturbations are decomposed into the adiabatic field Q a and the entropy field Q s 
as 

Q 1 = Q a e[ + Q s c\ ■ (13) 

Hereafter, for simplicity, we will concentrate on the two-field case although the extension to more fields is straightfor- 
ward. 

For the single field model, since the comoving curvature perturbation 1Z is given by 

^=(^=)^*> (14) 
the power spectrum of the primordial quantum fluctuation was given by [68j 

„ > , v 1 En 1 H 

p *' w =^77r^«- (15) 

where it should be evaluated at the time of horizon crossing c s „fc = a*iJ». The spectral index is 

d\nV n ,(k) 

n-n, -1= — -7j— t — = -2e-L-x- (16) 
However, as explained in the following sections, this will not be the case for multi-field models. 
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III. THE TRISPECTRUM FROM THE CONTACT INTERACTION 



In this and next sections, we will derive the four-point functions of the field perturbations at leading order in 
the slow- roll expansion and in the small sound speed limit. There are two important tree-level diagrams for the 
trispectrum. One is a diagram where the interaction occurs at a point, i.e. a contact interaction diagram and the 
other is a diagram where a scalar mode is exchanged. In the single field DBI inflation model, it was recently shown 
by [13, that both contributions are comparable. Because it is expected that this is also true for the two-field DBI 
inflation model, it is necessary to consider both contributions. In this section, we summarize the calculation of the 
four-point functions of the field perturbations at horizon crossing coming from the contact interaction diagram [64j . 
In the next section, we shall calculate the scalar exchange four-point function. 

Using the approximations mentioned above and following the ADM formalism [f| [T3, 0, [§8[ , the action up to fourth 
order can be calculated as [fijl 



S, 



(mc 
(2) 



1 



s. 



(main) 
(3) 



= - / ,//,/".,■ — 
1 

2 



Ql + Ql - ( diQ a d l Q a + d.Qsd'Q, 



= r / dtd 



Ql + Q a Q 2 s + -§ ( (diQ^Q, - d t Q a d l Q a ) Q a - 2 {diQ^Q,) Q t 



(17) 
(18) 



S, 



(4) 



— / <l.r''<ll- 
16 



c Ac 



2c 



5Ql + 6QiQi + Qt - ( SQidtQ^Qv - Q^Q s d l Q s + AQ a Q s d l Q cr d t Q s 



-Q%Qad l Qa + Q%Q s d l Q 



+^ (diQad l Q a ) 2 - 2 (d % Q«d l Q a ) (djQ.&Qs) + 4 (ftQ^Q.) + (diQs&Q,)' 



(19) 

It is worth noting that Eqs. (fT7|) . (fT8|) and (p~9|) can also be derived using a simpler method based on a Lorentz boost 
from the frame where the brane is at rest to the frame where the brane is moving 64 1. 

Based on these actions, we can obtain the third- and fourth-order interaction Hamiltonian densities as 16411 



njint 
"(3) 



n 



(4) 



2v/2AV 



4A c5 



2 

Ql + Q a Ql - % {Q ff (d t Q a d l Q, ~ d % Q s d l Q s ) + 2Q s ^Q CT 9 l Q s } 

2 

Qt + QlQl + J {{diQsd'Qs) Ql + {diQs&Qs) Ql) 



(20) 
(21) 



where here, Q a and Q s are the interaction picture fields. One important remark is that while the cubic part of H mt 
is the opposite sign of the cubic Lagrangian density C mt , this is generally not true at fourth order. 

From the second-order action, we can solve for the perturbations and quantize them according to the standard 
procedures of quantum field theory: 



^( r ?' x ) = ?i / [«n(l?.k)flnW +<(*7,-k)4(-k)] e lkx , 



(22) 



where a„(k) and 4( — k) are the annihilation and creation operator respectively, which satisfy the usual commutation 
relations: 

[a n (k!),c4(k 2 )] = (27r) 3 ( 5( 3 )(k 1 -k 2 )J nm , Mki), a m (k 2 )] = [4 (k x ), 4, 0*0] =0. (23) 

At leading order, solutions for the mode functions are given by 



u„(r], k) = An-^ (1 + *^ c s?7) e 
The two point correlation function is then obtained as 



—ikcst] 



{0\Q n (v = 0M)Qm(v = 0,k 2 )|0) = (2 7 r) 3 (5( 3 )(k 1 +k 2 )P Q „^(5„ m , 



(24) 



(25) 
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where the power spectrum is defined as 

^ = ^f> \A^ = \A S ? = ^N\ (26) 

and it should be evaluated at the time of the sound horizon crossing c s 

In terms of these quantum operators, the connected four-point correlation function coming from the contact inter- 
action diagram is given by [9^ | 

(fi|Q ro (0, k x )Q n (0, k 2 )Q p (0, k 3 )Q g (0, k 4 )|fi) c/ = 

*?<o| [Q m (o,k 1 )g n (o,k 2 )g p (o,k 3 )Q,(o,k 4 ), J ff ( ™ ) t (77)] |o) , (27) 

where Q m on the r.h.s. of the equation are the interaction picture fields and HY^f is the integral of the fourth-order 
Hamiltonian density ([21]), that is, = J cPxHjtf. 

The purely adiabatic, purely entropic and mixed components are respectively given by 



(rj|Q CT (o,k 1 )Q CT (o,k 2 )Q CT (o,k 3 )g g (o,k 4 )|r!) c/ = ^^(Vko J* 1 (-36^) 

i=l ZJi oc s J-J-i = i«i 



(28) 



n 4 a- 2 ^ 



H s 1/1 



(r!|Q s (o,k 1 )Q s (o,k 2 )g s (o,k 3 )g s (o,k 4 )|o) c/ = ( 27r ) 35(3) (E ki )^^nT^? {-g A ») > ( 30 ) 



''4 



. fc 2 fc 2 (k 3 -k 4 ) / 3(fc 3 + fc 4 ) , 12fc 3 fc 
= l 1 + K + 



perm. (31) 
H s -If,. 1 



(fi|g CT (o,k 1 )g CT (o,k 2 )g s (o,k 3 )g s (o,k 4 )|r!) c/ = (2^ (3 HE k i)^6nv( 6A i + 2 M ) > (32) 

(33) 



. _ fc 2 A: 2 (k 3 -k 4 ) / 3(fc 3 + fc 4 ) , 12fc 3 fc 4 

where "perm" denotes twent y th ree permutations of {fci, fe, fc 3 , fc 4 }. 

Following the analysis of [lOOj ] , if we describe the conversion of the entropy perturbation into the curvature per- 
turbation by a transfer coefficient T-rs, the final comoving curvature perturbation 1Z is expressed in terms of the 
adiabatic and entropy field perturbations as 



K = A a Q a * + A S Q S * , A a = (^7=) » A* = T ns j • (34) 

This implies that the final comoving curvature perturbation power spectrum is 

Vn = {l + Tl s )V n „ (35) 

where the horizon crossing power spectrum V-n, w & s given by Eq. (|15p . 

At leading order, the connected four-point function of 1Z coming from the contact interaction diagram is given by 1 

<^(k 1 )^(k 2 )^(k 3 )^(k 4 )) c/ = ^(Q CT (ki)g.(k 2 )g ff (k 3 )g ff (k 4 )) c/ 

+AlA 2 s ( (Q ff (ki )Q a (k 2 )Q s (k 3 )Q s (k 4 )) CI + 5 perm.) 

+^i 4 (g s (k 1 )g s (k 2 )g s (k 3 )g s (k 4 )) c/ , (36) 



There are other contributions to the connected four-point function of Ti, for example those coming from the five- and six-point functions of 
the scalar fields that can be expressed in terms of the power spectrum using Wick's theorem. In this work, we ignore these contributions 
because we are only interested in the intrinsically quantum four-point function. These extra terms are generated by the non-linear 
evolution outside the horizon and can be important 
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where "5 perm." denotes five permutations of the four- momenta. Using Eqs. (j28|) . (|30j) and (I32)) . it is possible to 
show a simple relation among the three different terms in the previous equation, as 

<g.(ki)Q CT (k 2 )g CT (k 3 )Q (T (k 4 )) c/ + <g s (k 1 )g s (k 2 )Q s (k 3 )g s (k 4 )) c ' / = (g ff (k 1 )g ff (k 2 )g s (k 3 )g s (k 4 )) c/ + 5 P erm. 

(37) 

This relation is perhaps less surprising if one notes that the different fourth-order interactions present in Eq. (j2Tj) 
also obey a similar relation, i.e. the purely adiabatic plus the purely entropic interactions (vertices) are equal to the 
mixed vertices. 

Using the previous results, the contact interaction four-point function of 1Z can be written as 

(^(k 1 )^(k 2 )^(k 3 )^(k 4 )) c/ = Ai(i + Tls) (^(g.(k 1 )g (T (k 2 )g ff (k 3 )g c7 (k 4 )) C7 

+Tl s (Q s (k 1 )Q s (k 2 )Q s (k 3 )Q s (k 4 )) CI ^j . (38) 

One immediately sees that if the transfer coefficient is non-zero, one can in principle use the trispectrum to distinguish 
single field from multi-field DB1 inflation. 



IV. THE TRISPECTRUM FROM THE SCALAR EXCHANGE INTERACTION 



In this section, we will calculate the four-point functions of the field perturbations at horizon crossing coming 
from a diagram where a scalar mode is exchanged. Within the "interaction picture" formalism [99j . to calculate the 
four-point function resulting from a correlation established via the exchange of a scalar mode, one needs to evaluate 
the following time integrals 



(fi|Q m (o, k 1 )Q„(o, k 2 )g p (o, k 3 )g 9 (o, k 4 )|0) 



SE 



di] 



di}(o\ g m (o,k 1 )g„(o,k 2 )g p (o,k 3 )g 9 (o,k 4 ),^ ( ™ ) t (r 7 ) ,H$(fj) |o) 



(39) 



where Q m on the r.h.s. of the equation are the interaction picture fields and -H| 3 * is the integral of the third-order 
Hamiltonian density flU, that is, Hjtt = J d 3 xH\™y 

The calculation of these time integrals is rather long and tedious. In Appendix^ we present an elegant and simple 
way of evaluating these contributions. This method is based on diagrams and simple rules that make the origin and 
physical meaning of all the different terms in the following expressions very clear. The purely adiabatic four-point 



function coming from the scalar exchange diagrams is identical to the single field one [62J , and it is given by 
(O|Q CT (0,k 1 )Q ff (0,k 2 )g ff (0,k 3 )Q CT (0,k 4 )|r!) 5£ = (2tt) 3 S^(K) u> 



16X c 9 s (kik 2 k 3 k 4 )i 
9( Fi(ki, k 2 , -k 12l k 3 , k 4l k 12 ) - Fi(-ki,-k 2 , -k 12 ,k 3 , k 4 ,k 12 ) 



3c s (k 3 • k 4 ) T 3 (k 4l k 2 ,-k 12 , fc 12 , k 3 ,k 4 ) - F${-k\, -k 2 ,-ki 2 , k 42 ,k 3 , fe 4 ) 



2(k 34 • k 4 ) T 3 {k 4 ,k 2 , -ki2,k a , fc 4 , k 12 ) - ^(-fci, -k 2 , -fcia, fc 3 , fc 4) k 12 



(ki • k 2 ) T 4 {-k 12 ,k u k 2l k 3 , fe 4 , fci 2 ) - :F 4 (-fc 12 , -k 2 ,k 3 , k 4 ,k 12 ] 



-2(ki2 • k 2 ) !F4,(ki,k 2 , -fa 2 , &3, foj, fcia) - ^(-fci, -&a, -fcia, &3, fc 4 , fei 2 ) 



(ki ■ k 2 )(k 3 • k 4 ) ( .F 2 (-fci 2 , fei, fc 2) k 42 ,k 3 , k 4 ) - ^(-fcia, -fa, -fc 2 , k 12 ,k 3 , k 4 ) 



-2(k x • k 2 )(k 34 ■ k 4 ) ^F 2 (-fe 12 , k x ,k 2 ,k 3 , k 4 ,k 12 ) - T 2 {-k 12l -k%, -k 2 , k 3 , k 4 ,k 12 ) 
-2(ki2 • k 2 )(k 3 • k 4 ) I T 2 (k 4l k 2l -k 12 ,k 12 ,k 3 , k 4 ) - J 7 2 (-h,-k 2 , ~ k 12 , k 12 , k 3 , k 4 ) 



+4(k i2 ■ k 2 )(k 34 • k 4 ) I T 2 {ki, k 2 , -k 42 ,k 3 , fe 4 , k i2 ) - T 2 {-k 4 , -k 2 , -k 12 , k 3 , fc 4 , k 12 ) 
-23 permutations of{fci, k 2 , k 3 , k 4 }, 



(40) 



where K = ki + k 2 + k 3 + k 4 and the information of the shape of the trispectrum is encoded in the T% functions 
defined by Eqs. (|B1|) . ()B2|) . (|B3|1 and (|B4[) in Appendix IB1 In the previous expression, "permutations" denote the 
other twenty three terms that result from the permutations of {fci, k 2 , k 3l k 4 } in the preceding term. 
The mixed component of the four-point function coming from the scalar exchange diagrams is 



(fi|Q CT (0, ki)Q ff (0, k 2 )Q s (0, k 3 )Q s (0, k 4 )|fi) S£ + 5 terms = (27r) 3 5 (3) (K) 



H 4 



lQX cl{kik 2 k 3 k 4 )2 



10 Ti{k\,k 2 , —ki 2 , fc 3 , k 4 , ki 2 ) - Fi(-kx,-k 2 , -fcia, k 3 , k 4 , k\ 2 



+2c 2 J - (k 3 ■^i){T 3 (k 1 ,k 2l -k 12 ,k 12 ,k 3l k 4 ) - T 3 {-k Xl -k 2 , -k 12 , k 12 , k 3 , k 4 



4(k 34 • k 4 ) T 3 {k 4l k 2 , -k 42 , k 3 , k 4 , k 42 ) - T 3 (-k 4 , -k 2l ~k 12 , k 3 , k 4 , fe X2 ) 



(ki • k 2 ) T 4 {-ki 2l k 4 ,k 2 ,k 3 , k 4 ,ki 2 ) - T 4 (~k 12l -k 4 , -k 2l k 3 , k 4 ,k 12 



+4(ki 2 • k 2 ) T 4 {k u k 2l -k 12l fc 3 , k 4 , k 12 ) - T 4 {-k 4l -k 2l -k 12 , k 3l k 4l k 12 



2c s (ki ■ k 2 )(k 3 • k 4 ) T 2 {-k 12 ,ki, k 2 , k 42 , fc 3 , k 4 ) - T 2 {-ki 2 , -ki, -k 2 , k 42 , k 3 , k 4 



+4(ki 2 • k 2 )(k 34 • k 4 )^F 2 (fci, k 2 ,~k 42 , k 3 , k 4 ,k 12 ) - T 2 {-k 4l -k 2 ,-k 42 , k 3 , fc 4 , k\ 2 ) 
+23 permutations of{fci, k 2 , k 3l k 4 }. 



(41) 



Finally, the purely entropic four-point function coming from the scalar exchange diagrams is 



(fi|Q s (0 ! k 1 )Q s (0,k 2 )Q s (0,k 3 )g s (0,k 4 )|r!) S£ = (2tt) 3 $< 3 >(K) : 



16X cf(fcifc 2 fc 3 fc4)s 
F\{ki, fa, -k 12 ,k 3 , k 4 , k 12 ) - Fii-ki, -fa, -fa 2 ,fa,fa, k 12 ) 



+c s (k 3 • k 4 ) ^{fa, fa, ~fa 2 ,ki 2 ,k 3 ,fa) - ^3(-fa,-fa, -k 12 ,k 12l k 3 ,k A ) 



2(k 34 • k 4 ) F 3 (ki,fa, ~ki2, fa, k t , k 12 ) - ^{-fa, -fc 2) -&12, ^3, fa, k 12 ) 



+ (ki • k 2 ) Ti{-ki 2 , k\,k 2 ,k 3 , fc 4 , fa 2 ) - Ti{-k\ 2 ,-kx, -fa, k 3 , fc 4 , fc i2 ' 



2(k 12 • k 2 ) T A {k x ,k 2 , -fa 2 , k 3 , k 4 , k 12 ) - ^{-fa, -fa, -fa 2 , fa, fc 4 , k 12 



c s (ki ■ k 2 )(k 3 • k 4 ) T 2 {-ki 2 , ki,k 2 , k\ 2 , k 3 , fa) - T 2 {-k 12 ,-k\, -k 2 , k\ 2 ,k 3 , fc 4 ) 



2(ki • k 2 )(k 34 • k 4 ) T 2 {-k 12 , fa,k 2 ,k 3 , fa,ki 2 ) - T 2 (-k 12 , -kx,-k 2 ,k 3 , fa,fa 2 ) 



2(k i2 • k 2 )(k 3 • k 4 ) T 2 {fa, fa, -fa 2 ,k 12 ,k 3 , fa) - T 2 (-fa,-k 2 , -ki 2 ,ki 2 ,k 3 ,fa) 



+4(ki2 • k 2 )(k 34 • k 4 )^ 2 (fci, fa, ~k 12 , fa, k A ,k X2 ) - T 2 (-k x , -fa, -k 12 , k 3l fa,k\ 2 ) 
-23 permutations of{fci, k 2 , k 3 , fa}. 



(42) 



The remaining possible four-point functions, (VL\Q <y Q r ,Q a Q s \Vl) SE and (Vt\Q (7 Q s Q s Q s \Vl) SE , are zero at leading order 
in the slow-roll expansion and in the small sound speed limit. 

As in the previous section, the connected four-point function of the curvature perturbation 1Z is related with the 
four-point function of the field perturbation as 2 

(^|^(k 1 )^(k 2 )^(k 3 )^(k 4 )|r!) S£ = ^<fi|Q CT (k 1 )Q CT (k 2 )g CT (k 3 )Q ff (k4)|0) S£; 

+A 2 a Al ((r!|Q CT (k 1 )g a (k 2 )Q s (k 3 )Q s (k 4 )|r!) S£ + 5perms) 

+^(o|Q s (k 1 )g s (k 2 )g s (k 3 )Q s (k 4 )|^) S£ , (43) 

Again, as in the case of the contact interaction trispectra, the scalar exchange trispectra obey 

(g CT (ki)g CT (k 2 )g CT (k 3 )g CT (k 4 )) S£ + (g s (k 1 )g s (k 2 )g s (k 3 )g s (k 4 )) SB = (g CT (k 1 )g CT (k 2 )g s (k 3 )g s (k 4 )) S£ + 5 P erm. 



(44) 



Using this fact Eq. (|43| can be simplified to 

(^(k!)^(k 2 )^(k 3 )^(k 4 )) SB = Ai(i + Tls) (^(g.(k 1 )g (T (k 2 )g ff (k 3 )g c7 (k 4 )) SB 

+T 2 <s (g s (k 1 )g s (k 2 )g s (k 3 )g s (k 4 )) S£ ^ . (45) 

Finally, the total four-point function of TZ, coming from the quantum four-point functions of the fields, is the sum 



2 Where again we ignore terms that can be expressed by copies of the power spectrum by using Wick's theorem. 
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of the contact interaction and scalar exchange contributions as 

(^(k 1 )^(k 2 )^(k 3 )^(k 4 ))* ota ' = (^(k 1 )^(k 2 )^(k 3 )^(k 4 )) c/ + (TZ(k 1 )n(k 2 )n(k 3 )TZ(k 4 )) SE 



K (1 + Tl s ) 



+T^ s f(Q s (k 1 )g s (k 2 )g s (k 3 )g s (k 4 )) c/ + (Q s (k 1 )g s (k 2 )Q s (k 3 )Q s (k 4 )) S£ 



(46) 



where the different terms can be found in Eqs. (|2"5)) . (jit})) . ([50)1 and (|4^|) . This constitutes one of the main results of 
this work. 

From Eq. (|46|) we can read the exact momentum dependence of the trispectrum which is related with the four-point 
function as 

<^(k 1 )^(k 2 )^(k 3 )^(k 4 )) totai = (2^) 3 ^ 3 )(K)T TC (k 1 ,k 2 ,k3,k 4 ). (47) 
In the next section we shall study the shape of the trispectrum in more detail. 

V. THE SHAPE OF THE TRISPECTRUM AND THE NON-LINEARITY PARAMETER t nl 

In this section, we will study the shape dependence of the different trispectra calculated in the previous sections. 
We will consider the so-called equilateral configuration. In the next subsection, we shall discuss some consistency 
relations that apply for the multi- field DBI trispectrum. After that, we will compute the non-linearity parameter tjvl 
and we shall plot the different contributions (pure adiabatic and pure entropy) to the total tnl both in the equilateral 
case and in a non-equilateral configuration. It is worth noting that from Eqs. (j5T|) and (j4"4"|) , the contribution from 
the mixed component can be given from the other two contributions. 

A. Consistency relations 

Here, we discuss whether the consistency relations obtained in the single field DBI inflation model can hold also 
for the multi-field DBI inflation model. In the following, as in 62], we will define the non-linearity parameter, which 
naively parameterizes the size of the trispectrum as 

4eC - ^~ ^t4 ,3 



(48) 



T A r L (ki,k 2 ,k3,k 4 ) = 1 H 2(i +T 2 ) ) TMki,k2,k 3 ,k4)II i=1 k 

(^1^2 + ^ffcf ) (^13 3 + ^ 14 ) + {k\k\ + fc 2 fcf ) (fc 12 3 + fc 13 3 ) + (&i A: 3 + fc 2 fc 3 ) (fc 12 3 + fc 14 3 



and we discuss behaviours of tnl- 

First, we discuss the consistency relation in the squeezed limit. As first pointed out by Maldacena Q and Seery et 
al. [3l| , all higher order correlators of 7Z in single field inflation obey consistency relations when one of the momentum 
vector is very small. Let us consider the limit ki — > 0, which implies that the corresponding mode 7?.(ki) leaves the 
horizon well before all other modes leave the horizon. By the time the remaining modes exit the horizon, 7?.(ki) will 
be frozen as a super-horizon mode and its only effect is to deform the background if the mode with ki which leaves 
the horizon earlier than the others is the adiabatic mode (<r). Then in this case, the trispectrum can be expressed as 

m 

(^(k 1 )^(k 2 )^(k 3 )^(k 4 )) -^(fci)-^-(^(k 2 )^(k 3 )^(k 4 )), (49) 

Hat 

where the power spectrum Vn(k) is given by (72.(ki)72-(k 2 )) = (27r) 3 <5( 3 )(ki +k.2)Vn(ki)- Using the fact that the bis- 
pectrum is Bq oc cJ 2 P^, it follows that in the squeezed limit the trispectrum is Tr oc (ecJ 2 P!^). Using P-jz oc (c s e) _1 , 
one finds that the total non-linearity parameter t^l should be at most of order ecj 2 (i.e. under our approximation it 
should vanish) when any of the momentum vectors goes to zero. From Eqs. (|28p. (150)) . (f52")l . (|4"0)l . (I4T]) and (|4"2")l . we can 
show that tjvl scales as 0(k 2 ), which means the consistency relation in the squeezed limit holds for the trispectrum 
in multi-field DBI inflation models. 
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However, if the mode with ki is the entropy mode (s), this explanation can not be applied. Instead, the fact that at 
leading order in slow-roll, the third-order action is written as Eq. (|18p and there are no three point functions like (sss) 
and (saa) turns out to be crucial. Actually, if one of the entropy modes leaves the horizon and becomes classical, the 
combination of the remaining three variables are (sss) for pure entropy component and (saa) for mixed component, 
respectively. As mentioned above, since these three-point correlations are absent at leading order in slow roll, the 
trispectrum vanishes in this limit. 

Similarly, we can discuss the consistency relation in the counter-collinear limit for the trispectra coming from the 
scalar exchange diagram. According to [38|, in the limit where the momentum of the scalar mode that is exchanged 
goes to zero, one can find a simple relation between the scalar exc hange tris pectrum and the power spectrum, in 
a similar way to Maldacena's consistency relations [f| (see also [28|, llOll Il02| |) since one can treat this mode as a 
background again if this mode is an adiabatic mode. Let us suppose that the momentum of the exchanged adiabatic 
particle is ki2 and that k±2 <C k\ ~ ^2,^3 ~ &4- Then the mode associated with this scalar particle will cross the 
horizon much before the other ki modes, where i = 1, ... ,4, and it only rescale the spatial background where the 
ki modes exist. Then, in the limit ki2 — > 0, the following relation should hold if the exchange mode is an adiabatic 
mode: 

<^(k 1 )^(k 2 )^(k 3 )^(k 4 )) SB -» (2ir) 3 S^(K)(n n - l) 2 Vn(k 12 )Vn(ki)V n (h) ■ (50) 

The previous equation implies that in the counter-collinear limit, the scalar exchange non-linearity parameter rjyf is 
of order e 2 (at leading order in our approximations this is equivalent to say that rjyf should vanish). From Eqs. (|40[) . 
(|4T|) and ([42]) . the non-linearity parameter t^l is shown to scale as 0(kf 2 ). Therefore, the consistency relation in 
the counter-collinear limit holds for the trispectrum coming from the scalar exchange diagram generated from the 
multi-field DBI inflation model. For the trispectrum coming from the other diagram, since there is no ki2 dependence, 
this consistency relation holds trivially. 

For the pure adiabatic, pure entropy and part of the mixed component, the exchanged mode is always an adiabatic 
mode. Thus we can apply the above argument to show that t^l vanishes in this limit. On the other hand, if the 
mode with ki2 is the entropy mode, which happens for the mixed component, we cannot apply the above arguments. 
However, in this case, the four-point function becomes a square of the two-point function of (sa). As was shown by 
Eq. there is no two-point correlation between the adiabatic and entropy mode at leading order in slow-roll. This 
explains the consistency relation in the counter-collinear limit for the remaining part of the mixed component. 

Now that we have checked the two important consistency relations, in the next subsection, we will concentrate on 
a concrete and best studied configuration in single field DBI inflation model, the so-called equilateral configuration 

Hums!. 

B. The shape of the trispectrum and tnl in the equilateral configuration 

In this subsection, we shall study the shape of the trispectrum in the limit where all the four momentum vectors have 
the same magnitude k. If we denote the angle between k^ and kj by 6ij, then in this configuration, we have cos(#i2) = 

C0S($34) = C0S($3), C0S($23) = COS^u) = COs(^i), COs(#i3) = COs($24) = COs(0 2 ) and COs(#i) + C0S(#2) + C0S(#3) = —1 

due to momentum conservation. Contrary to what happens in the case of the bispectrum, the equilateral configuration 
conditions do not fix all degrees of freedom (dof) required to describe the shape of the trispectrum and we are left 
with two angular dof (3l| . 

In the equilateral configuration, the contact interaction trispectrum is constant and only the scalar exchange 
trispectrum depends on the remaining angular dof, that we choose to be 9\ and 02- The plots of the shape of 
the scalar exchange trispectrum coming from the purely adiabatic component (the first term in Eq. (|43p) and the 
purely entropy component (the last term in Eq. (|4"5)l ) can be found in Fig. [TJ Using Eq. (|4"5)) . one can easily see that, 
for Tns <C 1 the purely adiabatic component dominates while for T-jis 3> 1, the purely entropy component dominates. 
Because the momentum dependence of the trispectrum changes depending on the value of T-jzs, the trispectrum can 
distinguish multi-field DBI inflation models from single field DBI inflation models if T-jzs is significantly large. The 
exact analytical expressions for the different contributions in the scalar exchange trispectrum can be found in Appendix 

D 

Even though the trispectrum includes all the information, for practical purposes it is more convenient to use another 
quantity t^l defined in Eq. (|48[) . In Fig. [2J we plot the non-linearity parameter r^|J for the purely adiabatic and 
purely entropy components calculated from the respective terms in the scalar exchange trispectrum, Eq. (|4"5)) . 

In Fig. [3 we plot the non-linearity parameter for the purely entropy component of the contact interaction 
trispectrum. The shape of the surfaces for the purely adiabatic and mixed components are the same as the one for 
the purely entropy component. However, the amplitudes are different, in particular, for the purely adiabatic plot, 
is negative. On the r.h.s. of Fig. [31 we plot the sum of the non-linearity parameters rj^£ and rjy^ that come 



11 




FIG. 1: The shape of the scalar exchange adiabatic (l.h.s panel) and pure entropy (r.h.s. panel) trispectra as functions of the 
variables cos Si and cos 62- The amplitude of the trispectra at the point cos 61 = cos #2 = — 1/3 has been normalized to unity. 




FIG. 2: Plots of the non-linearity parameter r^f calculated from the scalar exchange adiabatic (l.h.s. panel) and pure entropy 
(r.h.s. panel) trispectra as functions of the variables cos 9\ and cos 62- The amplitude of the first plot was rescaled by cf (1+T^ S ) 3 , 
while the amplitude of the r.h.s. plot was rescaled by c*(l + T^ s ) 3 /T^s ■ 

from purely entropic perturbations. If the transfer coefficient is large this will be the dominant shape in the total 
non-linearity parameter. 

Notice that since the factor kfj is multiplying the trispectrum, in terms of tnl , we confirm the consistency relation 
in the counter-collinear limit of the trispectrum coming from pure adiabatic, mixed and pure entropy components. 
When cos 6*i = — 1 or cos #2 = -1 or cos 83 = — 1 (on the diagonal lines in the plots of Fig. [21 when 0036*2 = — cos6*i) 
we can see the non-linearity parameters are zero. 

In order to quote a number for the value of the non-linearity parameter, we need to fix the remaining two degrees 
of freedom (#1 and #2) and in this way fully specify the momentum vectors' configuration. In the single field DBI 
model, there are good reasons to believe that the maximum of tnl occurs for the equilateral configuration with 
cos(#i) = cos(#2) = cos(#3) = —1/3 [62| 3 . In the multi-field case, for some contributions (like for example, Tjyf 
coming from the pure entropy trispectrum, as seen on the r.h.s. panel of Fig. [2]) the non-linearity parameter is not 
maximized for this particular angular configuration. However, for the sake of comparison with the single field model, 
we shall obtain the non-linearity parameter in that configuration. 



3 However, the definite proof that this is the case is still missing 
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FIG. 3: L.h.s. panel: Plot of the non-linearity parameter t nl calculated from the contact interaction pure entropy trispectrum 
as a function of cos#i and cos 62- R.h.s. panel: Plot of the sum of t^l and t^l coming from the pure entropy trispectrum as 
a function of cos Si and cos 62- The amplitude of the plots was rescaled by cf(l + T^ s ) z /T^ s . 



So, for the equilateral configuration, with cos(0i) = cos(6*2) = cos^) = —1/3, we obtain 

CI 0.036 CI 0.016T^s ci 0-052^ 

where the different values of tn l correspond to the three terms in the contact interaction trispectrum coming from 
the purely adiabatic, mixed and purely entropic components in Eq. (|36p . respectively. 
Similarly, from the scalar exchange trispectrum we obtain 

SE 0.60 SE o.67r|,g SE o.onr^g 

TNLi ° ~ ^(1 + 7^)3' TNL -ss ~ C 4( 1 + T 2 5)3 ' T NLis ~ 04(1 + 7^)3' (52j 

where again the different values of tnl correspond to the three terms in the scalar exchange interaction trispectrum 
coming from the purely adiabatic, mixed and purely entropic components in Eq. (|43p . respectively. 

It is worth mentioning that from Eqs. (f3"T)) and ([33]), both for the contact interaction (Eq. (|5Tj) ) and the scalar 
exchange interaction (Eq. (|52p). we can obtain the relation 

TNLaass ^ T^_ s T NLia + —^T NL4s . (53) 

Therefore, we conclude that the total non-linearity parameter tnl for the multi-field DBI inflation model in the 
equilateral configuration with equal angles between the momentum vectors is 

tot 0.56 0.687^ 0+27^ 

NLequi 4(1 + 2^)3 + 4(1 + 7^)3+4(1 + 7^)3- ^ ) 

It is worth noting that when there is no conversion of entropy perturbation into curvature perturbation, i.e. T^g = 0, 
Eq. reproduces the result obtained in the single field DBI inflation model [62}. On the other hand, when the 
transfer coefficient is large, i.e. Tug 3> 1, the final curvature perturbation is mainly of entropy origin and the last 
term in Eq. (|54p dominates over the others. 



C. The non-linearity parameter tnl in a non-equilateral configuration 

As we have mentioned before, the equilateral configuration is well investigated in the DBI inflation model. However, 
since there are five degrees of freedom to parameterize the shape of the trispectrum and the expression for the 
trispectrum is fairly complicated, non-equilateral configurations have attracted less attention. Because it has not 
been shown that the equilateral configuration maximizes the value of the non-linearity parameter, in this subsection 
we will consider a non-equilateral configuration. 

In our previous work [621 ] . we have studied the deviation from the equilateral configuration by fixing a particular 
set of degrees of freedom in the single field DBI inflation model case. Here, we apply the same approach and we study 
the behaviour of tnl in a non-equilateral configuration. 
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FIG. 4: Plots of the non-linearity parameter Tj^fJ calculated from the scalar exchange adiabatic (l.h.s. panel) and pure entropy 
(r.h.s. panel) trispectra for non-equilateral configurations as functions of the momentum amplitude q and cos 62- The amplitude 
of the first plot was rescaled by ci (1 + T^s) 3 , while the amplitude of the r.h.s. plot was rescaled by ci(l + T^ s ) 3 /T^ s . 




FIG. 5: Plots of the non-linearity parameter calculated from the contact interaction adiabatic (l.h.s. panel) and pure 
entropy (r.h.s. panel) trispectra for non-equilateral configurations as functions of the momentum amplitude q and cos 82. The 
amplitude of the first plot was rescaled by ci(l + T^ s ) 3 , while the amplitude of the r.h.s. plot was rescaled by Cg(l + T^ s ) 3 /T^ s . 

We consider configurations where the two momentum vectors IC3 and k4 have a different magnitude q from the 
magnitude of ki and k2 that we normalized to unity. Then in this configuration, we have cos(#34) = cos(#3), cos(#23) = 
cos(#i4) = cos(6'i), cos(#i3) = cos(#24) = cos(# 2 ), cos(#i 2 ) = q 2 (1 + cos(# 3 )) — 1. Momentum conservation implies 
cos #3 = — (1/q) (cos 9i + cos 6*2 ) — 1. We will allow one of the angular dof to vary freely, the other angular dof 6\ is 
fixed by cosB\ = —1/3 to the value that gives us the maximum for the single field model. In Fig. QJ we plot the 
non-linearity parameter rjyf calculated from the scalar exchange adiabatic (l.h.s. panel) and entropic (r.h.s. panel) 
trispectra, respectively. The maximum appears at q = 1 and cos(0 2 ) = —1/3, which supports the importance of the 
equilateral configuration. In Fig. [5J we plot the non-linearity parameter calculated from the contact interaction 
adiabatic (l.h.s. panel) and entropic (r.h.s. panel) trispectra, respectively. Comparing Figs. H and [5l one can see 
that around the maximum the purely adiabatic contribution coming from the contact interaction is smaller than the 
corresponding one coming from the scalar exchange interaction. Only for the entropy component, the amplitudes of 
the scalar exchange and contact interaction trispectra are comparable. We plot the sum of these two terms in Fig. [5] 
If the transfer coefficient is large this will be the dominant shape. 

VI. OBSERVATIONAL CONSTRAINTS 

Since we have obtained the full quantum trispectra for each component for the equilateral configuration with 
cos(#i) = cos(#2) = cos(6* 3 ) = —1/3, in terms of T-jis, we are in a position to discuss the relation between Tus an d 
observable quantities. It was shown in [53] that in the small sound speed limit {c 2 s ^ 1), the spectral index of the 
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FIG. 6: Plot of the sum of the non-linearity parameters and calculated from the pure entropy trispectra, for non- 
equilateral configurations, as functions of the momentum amplitude q and cos 82- The amplitude of the plot was rescaled by 
4{l + Tl s f/Ti s . 



scalar perturbations n-R in two-field DBI inflation model can be written as 



/3 l\j equil \ r f 

n n a -1— ^- -L- + a, sin(26) + 2/3* sin 2 9, (55) 

4 cos^ 8 Hj 

with a=-L, 0=*-z-^(lZ + %i), © = arcsin ( — ^= Y (56) 



aH' A 2 2 3H 2 y s a 2 )' \ y/T+T. 



where 1 and x are the slow- roll parameters defined in Eq. £ is the coupling between the adiabatic and entropy 
perturbation and fi s is the mass of the entropy perturbation, which are model dependent quantities. In order to 
obtain Eq. I|55[). they used the following relations: 

Vt 

r = — = 16ec s cos 2 e, (57) 
/^--i-cos^, (58) 



where r is the tensor to scalar ratio and /^f ™ parameterizes the size of the bispectrum in the equilateral configuration. 

In the ultra-violet (UV) DBI model (/ > 0) and without entropy transfer, i.e. Tns — (or equivalently sinG = 0), 
Eq. (|55p gives the well known stringent constraint as 



r> ; =(l-n n ), (59) 

fequi 
NL 



3.1|/^ U " 



by which one obtains r > 10~ 3 if one uses the observational constraints on /jy 9 ™ and 1 — n-ji from WMAP5 data 
[3(. On the other hand, since for general compact internal spaces, suggested by string theory, there is a geometrical 
constraint that limits r to be r < 10~ 7 the single field UV DBI inflation model is already excluded by the 
observational data. 

However, if we take into account the transfer from the entropy perturbation, the last two terms in Eq. (|55[) are no 
longer negligible. Furthermore, from Eqs. (|57p and (|58|. both r and /S^* are suppressed by cos 2 Q which scales as 
l/T^s f° r Tns ^> 1. Because of these factors, the constrai nt g iven by Eq. ([59]) cannot be applied and it is possible 
to satisfy the geometrical constraint r < 10~ 7 if Tns ^ 1 [53|- In such a case, since the last term of Eq. (fM)) will 
dominate over the other terms, t^™ can be expressed as r^?™ ~ 1.14(/^") 2 T^ S . Thus for a given fjq^ , t^vl ^ s 
enhanced for a large transfer T^s ^ U while the global magnitude of tjvl is reduced, going like 1/X^ 5 in this case. 



VII. CONCLUSIONS 



In this work, we have obtained the full quantum trispectrum of the primordial curvature perturbation at leading 
order in the slow-roll expansion and in the small sound speed limit in multi-field DBI inflation models. Previously, 
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we had obtained the contact interaction trispectrum [64|. In this work we completed the calculation of the quantum 
trispectrum including the scalar exchange trispectrum. 

We studied the momentum dependence of the trispectrum by separating it into the purely adiabatic, mixed and 
purely entropic components. For general configurations, as in the single field model, we confirmed that the consistency 
relations in the squeezed and in the counter-collinear limit hold. If a mode that crosses the horizon much before the 
other modes is an adiabatic one, these can be explained by the usual argument that the mode only changes the 
background (just like in the single field models). However, if the mode that crosses the horizon first is an entropic 
one, the consistency relations hold because of the absence of the two-point and three-point correlations among the 
modes remaining under the horizon scale. Essentially, the absence of these correlations is attributed to the particular 
properties of the DBI action discussed in [64] . 

In the equilateral configuration, we examined the momentum dependence of the scalar exchange trispectra as 
functions of the two remaining angular degrees of freedom. We found that the shape of the adiabatic, mixed, and 
purely entropic trispectra are different. This means that if the amount of the transfer from the entropy perturbation to 
the final curvature perturbation is significantly large, the trispectrum can distinguish multi-field DBI inflation models 
from single field DBI inflation models. In the equilateral configuration with cos 8\ = cos 82 = cos 83 = —1/3, we showed 
that for the contributions for tnl coming from the adiabatic and mixed trispectra, the scalar exchange diagram is 
more important than the contact interaction diagram, while for the pure entropy contribution both diagrams are 
equally important. This suggests that for quantitative discussions, the analysis done in Ref. (55| taking into account 
only the scalar exchange diagram is insufficient, at least, in this configuration. We have also considered configurations 
where the equilateral conditions are not satisfied. 

After obtaining the complete predictions for the quantum trispectrum in the multi-field DBI inflation model, we 
discussed the observational constraints on the non-linearity parameter, for the equilateral configuration mentioned 
above. In single field DBI inflation, some models based on string theory set-ups suffer from a severe constraint related 
with the geometrical structure of the compact internal spaces. However, in the presence of the transfer from the 
entropy mode, this constraint no longer holds. Furthermore, since both the tensor-to-scalar ratio and non-linearity 
parameter for the bispectrum are suppressed like T^J for a large transfer from the entropy mode to the curvature 
perturbation, it is easier to satisfy the corresponding constraints [5j| . On the other hand, the trispectrum is enhanced 
like i~]y 9 ™ oc T^tis^'nT 2 f° r a gi ven amplitude of the bispectrum ff* , while the global magnitude of tnl is reduced, 
going like l/T^ s in this case. 

In order to obtain a concrete value for the transfer coefficient Tpis , we need a specific model. It would be interesting 
to study this mixing in specific string theory motivated models [1031 ]. 

Note added: On the day this work appeared in the arXiv, the paper |l04| was also submitted to the arXiv, it 
also calculates the trispectrum in multi-field DBI-inflation. While we have concentrated on the trispectrum coming 



from the intrinsically quantum four-point functions of the fields which are generated during inflation, Ref. 1041 ] . by 
assuming a specific model, showed that the contribution to the trispectrum coming from the super-horizon non-linear 
dynamics is also important. 
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APPENDIX A: 



THE DIAGRAMMATIC APPROACH TO THE TRISPECTRUM FROM THE SCALAR 
EXCHANGE INTERACTION 



The calculation of the integrals in Eq. Q39p is rather long and it is easy to forget some terms or permutations. In this 
appendix, we present a systematic and simple way to keep track of all the different contributions and permutations 
to the integrals in (|39p . This approach is based on diagrams and simple rules that allows one to read the respective 
contributions for the scalar exchange four-point function directly from the diagrams. 

To calculate the scalar exchange four-point function we need to evaluate 



(fi|Q m (0, ki)Q„(0, k 2 )Q p (0, k 3 )Q g (0, k 4 )|fi) 



SE 



d V / dv(0\ Q m (0,k 1 )Q„(0,k 2 )Q p (0,k 3 )Q 9 (0,k 4 ),H ( - ) t (^) ,H$(fj) |0), (Al) 
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where on the l.h.s. of the equation Q m are the interaction picture fields and H%jf is the third-order interaction 
Hamiltonian that can be found in Eq. (|20|) or 

H$ (r?) = HT l (v) + H?* fa) + H™* (v) + H\ nt ( v ) + ( v ) , ( A2) 

where the different type of interactions are 

H{ nt ( V ) = J d 3 xg 1 aQ' rT \ W 2 nt (v) = J d 3 xg 2 aQ'„(dQ a ) 2 , H* 3 nt ( V ) = J d 3 xg 3 aQ' a Q' s \ 

Hl nt ( V ) = [ d 3 xg 4 aQ' a (dQ s ) 2 ,W 5 nt ( V )= [ d 3 xg 5 aQi s (dQ„)(dQ s ) 7 (A3) 



where here Q a and Q s are the interaction picture fields. The coupling constants are 

3 i = M, g 2 = -c 2 s M, 53 = M, g 4 = c 2 s M, g 5 = ~2c 2 M, (A4) 
where M = , 1 . The integrand of Eq. (| Al[) can be written as 

[[O™(0, k 1 )Q„(0, k 2 )Q p (0, k 3 )Q 9 (0, U),H$ (V)] , H$(rfi 

= J2 {FHi nt (rj)H} nt (fj) + H} nt (fj)Hi nt (r,)F - HJ 1 * (fj)F if*"* (rj) - Hj nt (r))F H* nt (fj) | (A5) 
i,j=i •> > 

where we have defined F as F = Q m (0, ki)Q n (0, k2)Q p (0, k3)Q g (0, k4). 

One can immediately see that at the leading order the only non-zero four-point functions 
are (Sl\Q*(0, ki)Q ff (0, k 2 )Q CT (0, k 3 )g CT (0, k 4 )\fl) SE , (O|Q s (0, ki)Q fl (0, k 2 )Q a (0, k 3 )g CT (0, k^)\fl) SE and 
(fi|Q s (0, ki)Q a (0, k 2 )Q s (0, k 3 )Q s (0, k A )\n) SE . The four-point functions {Q\Q a (o, k^Q^O, k 2 )Q a (0, k 3 )Q s (0, k 4 )\n) SE 
and (r2|Q<j(0, ki)Q s (0, k2)Q s (0, k 3 )Q s (0, k4)|il) s - E are zero at this order. This is a direct consequence of the absence 
of the interactions Q 3 and Q 2 Q S in the third order interaction Hamiltonian (|A2[) . In terms of diagrams this fact is 
clear. With the interactions (|A3|) one cannot draw diagrams with one and three external entropy legs, respectively. 

In order to find the non-zero four-point functions we need the following set of rules. We draw the allowed diagrams 
with two of the vertices (|A3|) , where one leg of the first vertex is glued with one leg of the other vertex (entropy legs 
can only be glued with entropy legs, similarly for adiabatic legs). The resulting four-point function can be constructed 
as: 

• Multiply by — 2(27r) 3 <!>( 3 )(K)iV 4 /(fcifc2fc3&:4) 3 / 2 , where N is the normalization of the mode functions. Multiply 
by the coupling constants of the vertices that enter the diagram. The overall minus sign is already present in 
(|Aip . The Dirac delta function ensures momentum conservation. The factor N 4 / 1 {kik 2 k 3 k4) 3 / 2 comes from the 
normalization of the four external legs, described by F . Finally, the factor of two can be understood if one notes 
that the contribution to the four-point function coming from the first term in the r.h.s. of Eq. (|A5p is real and 
the operator in the second term is just the Hermitian conjugate of the first term. Similarly with the third and 
fourth terms in the r.h.s. of Eq. (| A5|) . 



To each leg of the vertex with a spatial derivative associate the usual factor — ik. 

The different ways of drawing the same-looking diagram (ignoring the k labels but with the order in the left- 
to-right of the vertices the same) give a multiplicative symmetry factor. 

The contribution coming from the third (and fourth) term in the r.h.s. of Eq. (|A5p is equal to the contribution 
coming from the first (and second) term but with the sign of k\ and ki (but not fci2), that enter the arguments 
of Ti functions, changed to minus sign. 

In the diagrams with vertices H[ nt or H 3 nt the form factor is made of T\ functions. The definitions of the Ti 
functions can be found in Appendix [B] In the diagrams with vertices H 2 nt , Hl nt or H£ nt the form factor is 
made of Ti functions. In the diagrams where the l.h.s. vertex is either one of £/| n< or iJ™* and the r.h.s. vertex 
is either one of H 2 nt , H\ nt or H" Lt , the form factor is made of T3 functions. In the diagrams where the l.h.s. 
vertex is either one of H 2 nt , H" Lt or H™ 1 and the r.h.s. vertex is either one of H[ nt or H 3 nt , the form factor 
is made of T\ functions. The order of the six arguments in Ti is read from the diagrams too. The first three 
arguments are read from the l.h.s. vertex and the last three arguments from the r.h.s vertex. For example, for 
the diagram made of two iTJ vertices, the first and fourth arguments of Ti are the legs with time derivatives, 
ki in the first vertex and k4 in the second vertex. 
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FIG. 7: The two leading order vertices that contribute for the scalar exchange trispectrum in the single field case. Dot over 
a denotes that in that leg the interaction is a time derivative. V<r means that the interaction is a spatial derivative. There is 
momentum conservation at the vertices, hence at each vertex the momentum vectors have to sum to zero. Overall momentum 
conservation implies that ki2 = — k34, hence the labels in the r.h.s. diagram. 




FIG. 8: The two diagrams that contribute to the adiabatic four-point function, where the l.h.s. vertex is H[ nt and the r.h.s. 
vertex is H™*. The first, on the l.h.s. panel, we glue two legs with a time and a spatial derivative. The second, on the r.h.s. 
panel, we glue two legs with time derivatives. 



• At the end, add twenty three permutations of {k%, fe, ka, ki} to restore the symmetry of the four-point functions 
under the exchange of the momentum vectors. 

Let us consider an explicit example and calculate the contribution for the adiabatic four-point function of a diagram 
where the l.h.s. vertex is Iff 1 * and the r.h.s. vertex is H™*. These vertices are schematically depicted in Fig. [7] and 
they are the same ones as in the single field case. 

There are two distinct diagrams that one can draw, as can be seen in Fig. [5J hence there will be two terms. For 
the diagram on the l.h.s. of Fig. [8jwe get a symmetry factor six. This is because there are three ways we can glue the 
Vcr leg to the a leg and there are two Vct legs. For the diagram on the r.h.s of Fig. [8jwe get a symmetry factor three. 
This is because there are three ways of gluing the & leg of the H™ 1 vertex with one of the & legs of the H{ nt vertex. 
In both diagrams, we have two spatial derivatives so we get a factor (— i) 2 . These diagrams will be proportional to 

Finally, the diagram on the l.h.s. gives 



(n\Q a (0, ki)Q ff (0, k 2 )Q CT (0, k 3 )Q CT (0, k4)|fi) 



SE 



D (-«) 2 6(27r) 3 (5(K)- glg2 , (ki2 • k 4 ) {J 7 3(ki,k 2 , -k 12 ,k 3 , ki, k 12 ) - J 7 3{-k 1 ,-k 2 , -k 12 ,k 3 ,k 4 , k 12 )) 

[kik 2 k 3 ki) 2 

+23 permutations of{fci, k 2 , k 3 , k^}, (A6) 

and the diagram in the r.h.s. gives 

<r!|Q CT (0,k 1 )Q CT (0,k 2 )g ff (0,k 3 )Q ff (0,k 4 )|r!) 5£ 

D (-i) 2 3(27r) 3 c5(K)- ^TTs ( k 3 • k 4 ) {F 3 (ki,k 2 , —kn, k 12 ,k 3 , ki) - ^(-fci, -fc 2 , -fcia, k 12 ,k 3 , fc 4 )) 

(kik 2 k 3 k 4: )2 

+23 permutations of {ki, k 2 ,k 3 ,k4} (A7) 

Now, let us consider another example and calculate the contribution for the mixed four-point function of a diagram 
where the l.h.s. vertex is JJ|™* and the r.h.s. vertex is H^ nt . These diagrams involves vertices that contain entropy 
legs, they are schematically depicted in Fig. O 

Because the entropy and the adiabatic mode functions are the same, if one substitutes s with a in the vertices of 
Fig. [5] we get the vertices of Fig. \7\ Because the J-i functions only depend on the mode functions, that are the same 
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VS 



Vs 



Va 



FIG. 9: The leading order vertices containing entropy modes that contribute for the scalar exchange trispectrum in the multi- 
field case. Qi denotes the coupling constants. As before, the dot over s denotes that in that leg the interaction is a time 
derivative. Vs means that the interaction is a spatial derivative. A label s in the leg means that it is a entropy mode, a 
denotes an adiabatic leg. 




Vs s 





-k„ 



Vs Vs 




FIG. 10: The two diagrams that contribute to the mixed four-point function, where the l.h.s. vertex is HI and the r.h.s. 

jin 



for the entropy and the adiabatic modes we don't need to define other Ti functions, and the single field ones are 
enough. 

There are two diagrams that contribute in this example, as can be seen in Fig. 1 101 hence there will be two terms. 
In both diagrams, we have four spatial derivatives, this gives a factor of (— i) 4 . The diagram on the l.h.s. of the figure 
has a symmetry factor two, because there are two ways of gluing the s leg of the r.h.s. vertex with one of the two Vs 
legs of the l.h.s. vertex. Similarly, the symmetry factor of the r.h.s. diagram is also two. 

The l.h.s. diagram contributes with 

(O|Q CT (0, k!)Q CT (0, k 2 )Q s (0, k 3 )Q s (0, k 4 )|ft) S£ + 5 perms. 

D H) 4 2(2^) 3 ,5(K) ~ 2N g4g5 3 (k 2 • (-k 12 )) (k 3 • k 4 ) (^ 2 (fci, ha, -fci 2 , fci 2 , k 3 , k 4 ) - ^ 2 (-fci, -fc 2 , -fci 2 , fci 2 , k 3 , fc 4 )) 
(K1/C2K3K4) 2 

+23 permutations of{fci, k 2 , k 3 , fe 4 }, (A8) 
and the diagram in the r.h.s. gives 

(Sl\Q*(0, ki)Q CT (0, k 2 )Q s (0, k 3 )Q s (0, k 4 )\n) SE + 5 perms. 

D (-i) 4 2(2^) 3 ( 5(K)-^-^ r (k a • (-k 12 )) (k 12 • k 4 ) (JT 2 (fci, k 2 , -k 12) fc 4 , k 3 , k 12 ) - T 2 {-k u -k 2 ,-k 12 , k 4 , k 3 , k 12 )) 
{k\k 2 k 3 ki)2 

+23 permutations ofjfci, k 2 , k 3 , fct}, (A9) 

where 5 perms, denotes the following permutations of the displayed term [1 <-> 3], [1 <-> 4], [2 « 3], [2 <-> 4] and 
[!<-*• 3, 2 <->4]. 
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APPENDIX B: DEFINITIONS OF THE Ti FUNCTIONS 



In this Appendix, we present the definitions and the analytical expressions of the four Ti functions (with i = 1, . . . , 4) 
that appear in the main text. 

Because there are only four different types of double time integrals over the mode functions, we define four Ti 
functions (with i = 1, . . . ,4). Their explicit expressions are 



J — oo J — oo 



) 



6„6 



! 1 / A A 2 ' 



(Bl) 



where N is the normalisation of the mode functions given by N = H / \/2c s , A is defined by the sum of the last 
three arguments of the Ti functions as A = k 4 + k 5 + fc 6 and C is defined by the sum of all the arguments as 
C = ki + k 2 + k 3 + k 4 + k 5 + k^. The remaining functions are 



Mki,k 2 ,k 3 ,k 4 ,k 5 ,k 6 ) = [ d V a( v ) r dna{n)u*\nM)u%nM)u%nMU*\nM)u*{nM)u*{nM) 

J — oo J — oo 



N 6 c 2 Ifcifc, 



1 



H 2 \k 2 k 3 k 5 k 6 \^ AC 



^ k 5 + k 6 ^k 5 k e 



A 



A 2 



-\ (k 2 + k 3 + k 5 + k 6 + ^ ((k 2 + k 3 ) (k 5 + fce) + 2k 5 k 6 ) + 2 M^L±M 



+ q2 y k 5h + (k 2 + k 3 ) (k 5 + fc 6 ) + k 2 k 3 

+\ (k 2 k 3 (k 5 + k e ) + 2k 5 k e (k 2 + k 3 )) + 2 k2k ^ k& 

, 6 (i i a i 7 \ i 7 t /i . , x . k 2 k 3 k 5 k 6 \ k 2 k 3 k 5 k 6 
+ I fc 2 A)3 (fc 5 + fee) + fefce {k 2 + k 3 ) + 2 — I + 24 — 



(B2) 



T 3 (k\,k 2 ,k 3 ,k 4 ,k 5 ,k 6 ) 



o en , 

dijafa) / dfja(rj)U* (r), h)U* fa, k 2 )U* ( V , k 3 )U* (fj, k 4 )U*(fj, k 5 )U*(fj, fc 6 ) 

-oo J —oo 

N 6 4 \hk 2 k 3 k 4 \i 1 



H 2 \k,kS AC* 



fc 5 + ke. „fc 5 fc 6 3 /, , ^k^ke 
1 + . + 2^# + - fc 5 + fc 6 + 2 ■ 



A 



A 2 C 



A 



12- 



,fc 5 fc 6 



C 2 



(B3) 



T4,(ki,k 2 ,k 3 ,k 4 ,k 5 ,k 6 ) 



d V a( V ) / drjafajtf* fa, fci)*7*fa, Wfa, W fa, k 4 )U* (fj, k b )U* (fj, fc 6 ) 



= 2 



A 6 cf |fcifc 4 fc5fc 6 |^ 1 
H 2 \ k2 k 3 \l A*C 



A A 2 k 2 + k 3 



C C 2 



C 



+ 2 



A(k 2 + k 3 ) + k 2 k 3 
C 2 



A A 2 ' 
+3^ (A (k 2 + k 3 ) + 2k 2 k 3 ) + \2k 2 k 3 -^ 



(B4) 



where the function t/fa, k) is a modified mode function given by 



U(ri,k)=N— i T (l + tkc sV )e^ kc ^. 

\k\2 



(B5) 



If the sign of the argument k of U is positive then U is equal to the mode function, if the sign is negative then U 
equals the complex conjugate of the mode function. 
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APPENDIX C: THE SCALAR EXCHANGE TRISPECTRUM IN THE EQUILATERAL LIMIT 

In this Appendix, we present the scalar exchange trispectrum of the primordial curvature perturbation 1Z in the 
equilateral configuration. These analytical expressions can be derived from the general expression for the scalar 
exchange trispectrum, Eq. (|43[) . We will use these simpler equations to study the shape dependence of the trispectrum 
and to calculate the non-linearity parameter tnl in subsection IV Bl 

At leading order in the small sound speed and in the slow-roll expansion, the scalar exchange trispectrum is 



(fi|ft(0, ki)ft(0, k 2 )ft(0, k 3 )7e(0, k 4 )|0) SB = (n|ft(0, ki)ft(0, k 2 )7e(0, k 3 )ft(0, k 4 )|0) SB 

+(fi|7£(0, k x )7e(0, k a )ft(0, k 3 )ft(0, k 4 )|fi) 



in 

SE 

aass 



+(n\n(o,Kk 1 )K(OMMOM)K(o,k 4 )\n) SE 4s , (ci) 

where the three different contributions come from the terms in Eq. (|43|) that are coming from the purely adiabatic, 
mixed and purely entropy components. 
Explicitly, they are given by 

(n|ft(0, k x )7e(0, k 2 )ft(0, k 3 )ft(0, k 4 )|fi) 5£ 4ff = (27r) 3 ( 5< 3 > (K) 1 1 Vi (k 12 /k)V^ (k)V n * (*ia ) + 2 perms, 

c s U + i TC<sJ 

(C2) 

(O|ft(0,ki)ft(0, k 2 )ft(0,k 3 )7e(0,k 4 )|ft) SB CTass = (2^ 3 >(K) jfg , 3 4 p 2(fci 2 /fc)^(fc)^(fc 12 ) + 2perms, 

(C3) 

(n|7J(0,ki)%(0,k 2 )W(0,k 3 )7e(0 ) k4)|n) SB 4 , = (2ir) 3 6^(K) j|g , 3 ^P 3 (fci 2 /fc)^(fc)^(fc 12 ) + 2 P erms, 

(C4) 

where we have defined three T>i(x) functions, with i = 1, 2, 3, as 

x 4 (61952 + 48384a; - 10240a; 2 - 20832a; 3 - 5184a; 4 + 1412a; 5 + 824a; 6 + 103x 7 ) 



V 2 (x) 



2 9 (x + 2) 4 

c 4 (31232 + 23744a; - 6272a; 2 - 10496a; 3 - 1248a; 4 + 1904a; 5 + 824a; 6 + 103a; 7 ) 



2 8 (x + 2) 



4 



_ , . a; 4 (512 - 896a; - 2304a; 2 - 160a; 3 + 2688a; 4 + 2396a; 5 + 824a; 6 + 103a; 7 ) 

V3{x) = n^rw ' (C5) 



where fc 42 = ky/2(l + cos 6* 3 ), fci 3 = ky / 2(l + cos 62) and fci 4 = ky/2{l + cos Si), k denotes the common amplitude 
of all four momentum vectors kj, with j = 1, ... ,4 and the rescaled power spectrum is Vn*(k) = H 2 / (4ec s fc 3 ) that 
should be evaluated at sound horizon crossing. 

It is worth mentioning that from Eqs. (|3T[) and (|44|l . the following relation holds, 



V 2 =V 1 +V 3 . (C6) 
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